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ABSTRACT 


This report contains the development and utilization of a computer 
simulation of a planar phased array of circular waveguide elements. 

The simulation allows experimental work to be done on the computer and 
is directed toward the investigation of mutual coupling and wide angle 
impedance matching in phased arrays. Special emphasis is given to 
circular polarization. The aforementioned computer program has as 
variable inputs: frequency, polarization, grid geometry, element size, 

dielectric waveguide fill, dielectric plugs in the waveguide for 
impedance matching, and dielectric sheets covering the array surface 
for the purpose of wide angle impedance matching. Parameter combinations 
were found which produced reflection peaks interior to grating lobes, 
while dielectric cover sheets were successfully employed to extend the 
usable scan range of a phased array. The most exciting results came 
from the application of computer-aided optimization techniques to the 
design of this type of array 0 Different combinations of parameters 
were allowed to vary and several optimizations were found. As a conse- 
quence of the many separate investigations made possible by this program, 
some fairly general conclusions regarding polarization and WAIM effects 
could be made, and these are summarized. 

Also included in this report are a comparison of the major methods 
for wide angle impedance matching, a discussion of multimode elements, 
and a coverage of bandwidth limits set by not having true time delay 
phase shifters in an array. 
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I. INTRODUCTION 


As is well known, the variations of antenna element impedance 
with scan angle in electronically steered phased arrays can seriously 
affect radiation power efficiency with the problem being especially 
troublesome for wide angle scans. If an array antenna element is 
impedance matched to its source for broadside pointing, then any 
variations in its imepdance due to ever present mutual coupling 
effects as the beam is scanned will result in available power not 
being radiated, but reflected back toward the array generators. As 
a consequence, whenever there is a premium on power or where there 
is a need for wider efficient scan ranges, it is important that the 
array be optimally impedance matched for all scan positions. It 
was the objective of this work to investigate and expand upon methods 
for improving the efficiency of phased array antennas. Most of the 
compensation measures considered work by providing a type of feedback 
which continuously adjusts so as to offset element impedance and conse- 
quent efficiency changes with scan angle. At the time of this study 
most of the existing work was concerned with linearly polarized sources, 
and so prime effort was given to extending the compensation procedures 
to matching of circularly polarized elements. The results of this 
endeavor are contained in three technical reports and this final 
report. 
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Technical Report No® 1, "Antenna Arrays: Efficiency and Effi- 

ciency Improvement Through Compensation," [1] is a presentation and 
analysis of the major compensation techniques for wide angle impedance 
matching (WAIM) ® The procedures described are generally separated into 
one of two categories ® The first considers methods which rely on 
structures which modify the array surface environment® Examples of 
such are baffles and dielectric cover sheets [2]® The other category 
contains those methods which utilize matching schemes located behind 
the array surface, examples of which are Hannan’s interconnecting 
circuit method [3] and the multimode element approach [4]. The multi- 
mode element procedure uses internal irises and discontinuities to 
generate higher order element modes which can hopefully be adjusted to 
offset mutual coupling scan variations® Since this method was omitted 
from Technical Report No® 1, it is included in Chapter II of this 
report® Also in Technical Report No® 1 is a good survey of design 
results for linear and planar phased arrays which should be useful to 
anyone needing a single-source compilation of phasing, grating lobe, 
and grid geometry considerations in phased arrays® A thorough coverage 
of array efficiency is given with two levels of approach® One treats 
all elements as ideal and determines maximum possible array efficiency 
as a function of grid geometry and spacing, and the other considers 
mutual coupling effects and illustrates the steps involved and the 
problems associated with the analysis of arrays of physical elements® 

Technical Report No® 2, "Microstrip Technology and its Application 
to Phased Array Compensation," [5] started out as a specific investi- 
gation of the practicality of Hannan’s interconnecting circuit method 
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for wide angle impedance matching and led to the need for a fairly 
thorough review of microstrip fabrication and design technology, A 
summary of these findings occupies the first portion of that report* 
Materials, element realizations using "lumped" versus "distributed" 
methods, processing, expected frequency and cost limits, losses, and a 
determination of microstrip characteristic impedance which better 
matches experimental measurements for lines with small o>/d ratios are 
all included in early sections. With the preceeding as a basis, cal- 
culations are made for wide angle impedance matching an array of 
S-band, half-wave dipoles using Hannan’s method* The design is carried 
all the way to the point of' specifying size and component construction 
for the interconnecting networks using MIC techniques. It was found 
that the required compensating networks were realizable for the 
S-band case using lumped MIC deposited elements* In order to obtain 
some of the relatively small effective inductances and capacitances 
required for matching, off-resonance, parallel LC tank circuits were 
used so that larger value.* better quality inductors and capacitors 
could be constructed* 

Based on our study the following observations are made relative 
to the usefulness of the interconnecting circuit method of compen- 
sation* 

Disadvantages : 

1. This method requires "a priori" knowledge of array impedance 

variations with scan angle usually necessitating array simulation 
or experimental measurements. The later can be rather expensive 
if an array is built and then it is found that it cannot be 


matched satisf actorily « 
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2o Presently this technique is not cost competitive with seme of 
the ether matching methods* especially . in terms of the cost of 
design time. Often a few crude calculations and some simple 
measurements can achieve a reasonable wide angle scan improvement 
using the dielectric cover sheet method,, whereas great care and 
much time must be devoted to obtain a possibly temperamental 
microstrip realization:, The real hope for application of Hannan’s 
procedure would be in the future for phased arrays of the modular 
type which have, active sources behind each antenna element. If 
subarxays of many such elements could be batch processed, then 
the production cost with interconnecting -circuits would be 
essentially the same as without them 0 It would still be necessary 
to know mutual coupling and impedance effects prior to the subarray 
synthesis so that correct compensating circuitry could be specified. 
Although it is possible to individually trim values of L and C 
after the elements are fabricated, it would be foolish to rely on 
such adjustments for arrays with hundreds of elements 0 Instead 
it is desirable to take adequate steps to make accurate predictions 
of needs. Note the advantage of having computer-aided simulations 
describing array performance. 

3. Matching networks might limit the frequency response of arrays 

intended for wideband operation. Five percent bandwidths, however, 
should not show significant degradation G 

Advantages : 

I. Extremely good matches are theoretically possible [6l> in fact 
much .better matches could probably be obtained than with any of 


the other WAIM schemes. 
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2 0 With individually adjustable compensating circuits, it would be 
possible to make pos t-f abrication system alterations. Also the 
capability of changing network values in the neighborhood of the 
array edge might alleviate some potential edge-effect problems. 
Finally, Technical Report No 0 2 contains many of the fabrication 
related problems and experiences relating to other facets of the 
contract program^ For instance, a circularly polarized, Ku-band 
array of planar spirals was photcli thcgraphically produced with 5 mil 
wide gold conductors on an alumina substrate 0 Broadband baluns for 
feeding wideband elements were investigated, built, and res ted 0 In 
addition, power limits set by using microstrip techniques and materials 
to realize physically small, Ku-band spiral antennas were estimated „ 

The spiral antenna element is analyzed in Technical Report No 0 3, 
"Simulation of Planar Spiral Antenna Elements for Phased Arrays c " With 
reliable computer-aided simulations and designs valid for large arrays 
of spiral antenna elements, future array systems using spiral elements 
could be made mere economical and realizable 0 The importance of this 
will be especially great for operation at X-band and above where elements 
are physically very small (for half-wavelength spacings between elements 
at 15 GHz the spiral diameter has to be less than loO cm), and as a 
consequence, . exacting. tolerances are required in construction,, To 
meet these standards all the elements in a high frequency array could 
be formed in cne operation using solid-state and microstrip deposition 
techniques o It seems reasonable to think that an entire "integrated" 
system including RF electronics, antenna elements, interconnecting 
circuits for wide angle compensation, phase shifters, and control cir- 
cuitry could be built as a single unit using this technology,, In any 
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event, systems requiring the above processes to produce the antenna 
elements are not recommended for experimental optimization since the 
resultant arrays are haphazard, expensive, and inflexible 0 Thus, if 
a sufficiently general computer program for simulating arrays of spiral 
antennas was part of the public domain, all necessary experimentation 
and optimization' could be computerized and the availability of active 
arrays for X-band and Ku-band might be considerably advancedc 

In Technical Report No, 3a computer-aided solution for mutual 
coupling and impedance effects in a spiral antenna is obtained, and 
groundwork is laid for an analytic solution for arrays of Archimedian 
spiral radiatorsc Solutions use Harrington’s method of moments [7] 
and a computer .matrix inversion. It was found that conductor dimensions 
in the spirals are more critical in- determining impedance than 
originally anticipatedo Further, it was noted that the reactive portion 
is not always negligible even when dimensions are nearly those to give 
complementary . symmetry o Some reasonably good designs are obtainedo 
Topics not covered in ether reports are covered in this final 
report o In Chapter II are two important phased array design consid- 
erations which relate to the problems of wide angle impedance matching. 
One is the use of multimode elements to achieve wide angle scan compen- 
sationo The other deals with bandwidth limits posed by not having 
true time delay phase shifters in a phased array 0 It is found for 
scans out to + 60 c on allowing a 0o7 db gain loss that the array beam- 
width in degrees is approximately equal tc the array bandwidth (%) for 
cw operation o 
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Chapter. Ill contains che development and utilisation of a 
computer simulation of a ..planar , phased ..array of circular waveguide 
elements o The simulation allows experimental work to be done on 
the computer and is directed toward the investigation of mutual 
coupling and wide-angle . impedance marching in phased arrays „ Special 
emphasis is given to circular polarization The aforementioned com- 
puter program has as variable inputs: frequency, polarization, grid 

geometry, element size, dielectric waveguide fill, dielectric plugs 
in the waveguide for impedance matching, and dielectric sheets covering 
the array surface for the purpose of wide angle impedance matching. 
Parameter combinations were found which produced reflection peaks 
interior to grating lobes, while dielectric cover sheets were success- 
fully employed to extend the usable scan range of a phased array . 

The most exciting results came from the application of computer-aided 
optimization techniques to the design of this type of array . Different 
combinations of parameters were allowed to vary and several optimi- 
zations were found- As a consequence of the many separate investigations 
made possible by this program, some fairly general conclusions regarding 
polarization and WAIM effects could be made, and these are summarized 
in Chapter III. 

The last chapter gives a brief description of some of the experi- 
mental work done during the course of this study. Of importance to 
the above was the design and construction of a waveguide simulator [8] 
which was used to test the accuracy of the computer simulation for 
the array of circular waveguides . 


CHAPTER II 


SUPPLEMENTARY PHASED ARRAY DESIGN CONSIDERATIONS 


In this chapter some topics are covered which are not included in 
other reports, yet are pertinent to wide angle impedance matching (WAIM) 
considerations in phased arrays. A brief survey and comparison of some 
of the major WAIM methods is presented below and later sections in 
this chapter are devoted to discussions of multimode elements for 
scan matching and bandwidth limitations in phased arrays. 

A. Methods for Wide Angle Impedance Matching - Comparison 

Two of the main objectives of this project were to determine design 
procedures for the various WAIM methods and to extend WAIM techniques to 
the matching problems associated with circularly polarized antenna elements. 

A thorough coverage of the theory and design practices for most of the 
WAIM methods is found' in Technical Report No. 1 [ 1] while Chapter III of 
this report contains an exhaustive study of wide angle scan improvement methods 
with circularly polarized elements. The analysis in Chapter III is for 
array elements described by circular apertures, cylindrical cavities, 
or circular waveguides. A treatment of the simpler analytic case of 
rectangular waveguides used for circular polarization is given in a 
recent article by Lee [9]. Technical Report No. 3 contains an investi- 
gation of Archimedian spirals. The above three works well represent the 
types of flush mounted elements commonly used for circularly polarized 
array elements. 
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Briefly the principal methods for wide angle impedance matching 

are : 

lo Magill and Wheeler dielectric cover sheet [2 ]: Here a 

dielectric slab with a high permittivity is placed over an array surface, 
and a height, thickness, and dielectric constant combination is sought 
which will provide WAIM. This method is inexpensive, lends itself 
nicely to all forms of investigation (waveguide simulators, computer 
simulation, and experimental coupling or impedance measurements), and as 
is shown in Chapter III is useful for matching arrays with either linearly 
or circularly polarized elements. Its only drawbacks are that matches 
usually are not ideal and bad parameter choices can worsen match rather 
than improve it. 

2. Edelberg and Oliner fences and baffles [11]: 

Metal fences or baffles or posts are placed on an array surface 
in order to reduce deleterious mutual coupling effects. In practice 
metal fences were used to reduce E plane scan variations in an array 
of A/2 dipoles. An array with a VSWR of <_ 2.9 out to a 60° scan 
angle had, after adding such fences , a VSWR <_ 1.6 [10]. In general 

this method relies on experimental trial and error. Also the surface 
modifications are not flush mounted and can contribute to shadowing 
which limits the potential scan range for the array. 

3. Interconnecting circuit method of Hannan, Lerner, and Knit tel 

[3 ]: 

This method places reactive networks between adjacent element feed 
lines in an array. The interconnecting reactances produce an equivalent 
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reactance to ground which varies with scan angle and can be used to 
counteract variations in element impedance with scan caused by mutual 
couplingo Technical Report No, 2 [ 5 ] contains an extensive discussion 
of this methodo Its main problem is that it is very difficult to 
implement, 

4o Multimode elements: 

These are discussed later in this chapter. If the weight of 
waveguide elements can be tolerated, this type of element is a 
good approach for scans out to 60°. They can be adjusted to impedance 
match the feed lines, are fairly broadband, and incorporate the means 
for WAIM into the element itself. They operate by generating higher 
order propagation waveguide modes which can be adjusted to cancel out 
reflections from the aperture. The' mutual coupling between elements 
for the various higher order modes tends to average out to give a 
wide matched scan. The concept might be usefully applied to other 
element types, 

5, Use of many closely spaced elements: 

This type of system can result in wide angle impedance match, 
however, it suffers' from two disadvantages . First the additional 
number of elements needed to fill the array add unnecessary cost 
compared to other matching schemes that can achieve the same match 
using the minimal number of elements. Finally very small cavity or 
waveguide-like elements are difficult to match to free space. 
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6. Electronically tunable matching networks: This is a futuristic 

approach which is presently expensive and subject to failure when components 
fail. 


B. Multimode Array Elements for Wide Angle Scans 

Wong et al. [4 ] have by generating higher order modes in elements 
for phased arrays achieved a novel method for controlling mutual 
coupling so as to maintain a nearly uniform driving point impedance 
over a specified scan range. Using rectangular waveguide elements 
enough propagating modes are generated (usually TE^q, ^20’ anc ^ 
and then adjusted to produce a cancellation in the net reflected 
field from the radiating aperture. The higher order modes are 
generated in waveguide inserts using discontinuities. These higher 
order modes accomplish two purposes. First the inserts can be adjusted 
to match the waveguide to the aperture radiating into free space (a 
broadside match is normally sought) . The second is a consequent 
side-effect which often improves wide angle scan performance. 
Qualitatively this results from the fact that the net mutual coupling 
impedance between a reference element and all the other elements in 
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the array now depends on several different modes whose changes with scan 
angle tend to offset one another, thereby maintaining a constant 
over wider scans. 

Several methods for generating the higher order waveguide modes 
are illustrated in Fig. 2-1. The two principal generating schemes 
suggested are discussed below. 

1. Guides are filled with a dielectric slab and a discontinuity 
is positioned to produce TE^q and TE^ modes in addition to the TE^q 
mode already present. A cylindrical hole is useful for such mode 
excitation. The length of the slab and the location of the holes are 
varied to obtain the relative magnitudes and phases of each mode 
necessary for matching. Waveguide simulatiors and computer simulations 
were used for initial designs, and at Hughes Aircraft [ 4] a 256 element, 
linearly polarized array of rectangular waveguide elements was built. 

The Hughes array had an input VSWR <_ 1.5 for scan out to ±47° in E, 

H, and D planes,- The VSWR was less than 2.5 for 60° scans. The 




Fig. 2-1. Multimode Waveguide Elements Realization. 
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major problem with the hole-in-dielectric type of element is one of 
bandwidth. Since the relative hole position changes with frequency, 
the mode balance is upset and the element is frequency sensitive. 

2. Inductive rises at the interface between waveguides and free 
space are capable of producing a wideband, multimode element. The size 
of an iris is found to be fairly effective in controlling the relative 
proportions of the higher order modes for match. In addition, no 
critical distances are invloved and thus the element has good frequency 
response characteristics. Arrays of such elements with a 20% bandwidth 
out to ±53° with a VSWR <_ 1.7 have been built and tested [ ^ ]. 

C. Bandwidth Considerations in Phased Arrays 

Whether or not bandwidth is part of the design problem in /phased 
arrays depends on the system objectives and on the frequency ranges 
used. For instance a 1% bandwidth about a 18 GHz carrier gives a 
180 MHz channel while at 3 GHz the cahnnel would be limited to 30 MHz 
(which is still adequate for many applications). The factors in- 
fluencing bandwidth are array aperture effects, the array feed system, 
the type of phase shifters, and component or receiver bandwidths. 

From a noise point of view it is often advantageous to let the 

antenna do the prefiltering for a receiver Irv which case narrow 

array bandwidths could be a desired asset. Emphasis in this section 

is on bandwidth limitations related to array apertures and feeds for arrays 

which are characterized by equal feed line lengths and fixed phase 

shifters. As a consequence of the fixed phase shifters, the array 

beam scans with frequency causing a loss in gain relative to the gain 
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expected for a given scan direction. Note that an array using true 
time delay phase shifters would not scan with frequency and would have 
an infinite array bandwidth. The bandwidth of true time delay system 
would be set by the receiver or individual components. 

For the fixed phase shifter array Frank and Cheston [29] arrive 
at the following easily remembered approximation for array bandwidth: 

1. For CW operation, a 60° maximum scan, and a 0.7 db allowed gain 
degradation 

r\j 

Bandwidth (%) = Broadside beamwidth (degrees) 

2. For pulse applications 

r\j 

Bandwidth (%) = 2x Broadside beamwidth (degrees) 


Two examples of this type of feed and phase shifter arrangement, the 
parallel feed and a modified corporate feed are illustrated in Fig. 2-2. 
The phase shifter in the corporate feed have been placed so that 0-360° 
phase shifters can be used without needing to drop multiples of 360° 



phase shifters 



Fig. 2-2. Array Feeris 
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For linear arrays such as shown in Fig. 2-2 the incremental phase 
shift between adjacent elements needed to steer the beam in the 0^ direction is 


27rf 1 d 

6 = kd sin0~ = sin0^ 

0 c 0 


( 2 - 1 ) 


where 


^ = kd sin0 - 6 and f is the center frequency. 


Providing the phase shifter settings remain constant, then at a new frequency 
f^ the beam is steered to the direction 0^ + A0 giving 

2irf d 27Tf d 

6 = — sin0 n = — sin(0 n + A0) (2-2) 

c 0 c 0 

For A0 small, sin (0^ f A0) ^ sin0^ + A0 cos0q. Substituting in the above 
2 tt f d a, 27Tdf 

1 sin0 Q * 2 (sin0 Q + A0 cos0 Q ) (2-3) 

c c 

and solving for A0 yields 

* f l- f 2 

A0 — — tan0Q (radians) . (2-4) 

From the above it is seen that for f^ greater than f^ the beam is scanned 
toward broadside, while for f ^ less than f^ the beam swings away from 
broadside. Further, when the beam is pointed toward broadside tan0^ = 0, 
and hence the beam does not scan with frequency (bandwidth infinite) . 

Additional simplifications result by assuming a small percentage 
bandwidth. Writing 


f 


2 


f i + 


Af 


(2-5) 
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it follows that 


f 1 -(f 1 + Af) _ Af tan0 

+ A0 = 1 _ tan0 n = + 0 

£ 1 + A£ £ 1 


( 2 - 6 ) 


or 

A6 (radians) = tane Q 

where 

BW = 2Af = bandwidth 
A0 = scan deviation from band center. 
In degrees Eq. 2-7 becomes 


(2-7) 


A0 (degrees) = .287 BW(%) tanO 


0 


( 2 - 8 ) 


Frank and Cheston define a bandwidth factor 


K = bandwidth factor = 


and normalize (2-8) to 


bandwidth(%) 


broadside beamwidth (degrees) (2-9) 


A0 ^ 

T— f rr - 0.287 K sin0 n 

0g (scanned) 0 

by introducing the large array approximation 


( 2 - 10 ) 


0 (scanned) = beamwidth at 0 n = 1 d - e ^ 

B 0 cos0q 


( 2 - 11 ) 


For a maximum gain loss of 0.7 db the bandwidth must be such that the 
beam never scans more then + 1/4 the local beamwidth and hence 


A0 


0 (scanned) 
n 


< 1/4 


( 2 - 12 ) 
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Applying this condition to Eq. 2-10 gives 

k = °4 

max sm0Q 

For an array designed of a maximum scan of 60° 


(2-13) 


K 

max 


.87 ^ , 

->o i 
sin 60 


Thus the worst case for a 60° from broadside scan is 


BW(%) = 0 (broadside) = broadside beamwidth (degrees) 

(2-14) 

A similar procedure working in terms of ^equivalent pulse length" 
results in 

O' 

pulse length ~ 2x aperture size 


bandwidth (%) = 2x broadside beamwidth (degrees) 


and 


CHAPTER III 


CIRCULAR WAVEGUIDE ARRAYS 


The main concern of this report is with antenna elements capable of 
producing circular polarization in a phased array environment and 
specifically with the mutual coupling and driving point impedance varia- 
tions with scan angle which can lead to severe power and array efficiency 
losses as the array is scanned away from broadside. Methods exist for 
wide angle impedance matching, but most of them require quantitative 
knowledge of the behavior of array coupling with scan in order to be 
applied. In this chapter analytical results and an effective digital 
computer simulation are developed which are valid for an infinite planar 
array of circular waveguide antennas, an important class of radiating 
elements useful for circular polarization. Using this simulation WAIM 
(wide angle impedance matching) methods are employed and computer-aided 
optimum array designs are determined using Rosenbrock's pattern search 
method [17J. In addition, these computer phased-array simulations will 
be used to hopefully resolve the controversy over two seemingly contra- 
dcitory ways to best achieve WAIM, namely, by using closely spaced large 
elements with a high degree of mutual coupling or by using small loosely 
coupled elements. 

The following antenna realizations are fairly representative of those 
used to obtain circular polarization in arrays: 

1. Circular waveguides or cavities with or without dielectric 
covers . 

2. Square guides. 
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3. Crossed dipoles (wire, thick, folded, bent)* 

4« Crossed slots backed with rectangular or circular cavities. 

5. Spirals. 

6. Wheeler elements [18] each consisting of a circular cavity below 
cutoff, fed with an internal crossed thick dipole, and having 
dielectric layers for transmission matching. 

Much theoretical and experimental work has already been done on rectangu- 
lar guide and cavity elements [36 -40] and on external dipole arrays [11-16], 
However, square and rectangular elements do not pack well in 60° -tri- 
angular grid arrays which are in standard use because of the savings in 
number of elements [21] , and further the element shapes are such that 
mutual coupling is apt to be greater in certain regions of an array cell 
than that found for circular elements. Volumes of data on dipoles are 
available, yet because dipoles are subject to shadowing with scan and 
because they require support structures which can interact with the dipole 
fields, dipole arrays can be more difficult tc analytically predict and to 
achieve wide angle scans. On the other hand, flush mounted circular ele- 
ments have the following advantages for use in phased arrays [18] : 

1. Circular apertures fit into equilateral triangular grids thus 
requiring the fewest number of elements for wide-angle scans 
free of grating lobes. 

2. Circular apertures are suited for resonant windows having 
maximum bandwidth for radiation loading, 

3. A circular aperture is resonant in the TE^ mode, 

4. Circular symmetry is suited for linear, crossed- linear , 
or circular polarizations. 

5. Circular symmetry is amenable to realization with waveguide 
simulators for infinite planar arrays, 

6. Circular apertures tend to reduce coupling wave effects, 

7. Mechanical strength and simplicity are characteristic, 

8. Many circular elements lend themselves to analysis. 
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9c Flush elements are aer ©dynamically smooth. 

In view of the above, it was decided that practical mutual coupling and 
efficiency effects in scanned arrays could be studied fruitfully by in- 
vestigating infinite planar arrays of circular waveguide elements and 
infinite arrays of planar spirals . Arrays of circular waveguides are 
treated in this report, while arrays of spirals are examined in Tech- 
nical Report Noo 3 [ 19]. A diagram of a possible realization for a 
typical radiator with circular geometry is shown in FIgr 3-1. 



Fig 


3-1 


Circular Radiator 
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As can be seen, the element in Fig. 3-1 consists of a circular waveguide 
with provisions for a feed and dielectric discs (some noncircular) for 
mode generation and impedance control* The noncircular discs can be used 
to generate and phase the modes needed to produce an arbitrary polarization, 
and the circular dielectric discs or "plugs'' can be used as resonators to 
obtain a wideband element or as quarter-wave plates for impedance matching. 
In any event, circular waveguide systems can be used to investigate a wide 
class of flush-mounted, circular aperture type sources (the circular cavity 
and the Wheeler element both easily fit this category) . The spiral is a 
class by itself, and the analysis of both types should further the state 
of the art and give insight into the optimum design of phased arrays. 

A. Computer Simulation of an Infinite Array of Circular Waveguides 

An infinite planar array of circular waveguides such as shown in 
Fig. 3-2 has been analyzed and simulated on a digital computer with the 
following design flexibility: 

1. Arbitrary parellelogram grids (a,b,d). 

2. Arbitrary dielectric constant in the waveguides (EPS). 

3. Dielectric plugs or discs (EPS2 , T2) in the waveguides for 
impedance matching. 

4o A dielectric cover sheet (EPS3, T3) over the array face for 
wide angle impedance matching (WAIM) . 

5. Scan in any direction (THETA, PHI) or across any scan plane* 

6o Arbitrary polarization using crossed TE^ waveguide modes* 

To solve this complex problem, Harrington's method of moments [7 ] 

or, specifically, the Ritz-Galerkin method is applied to a boundary value 
formulation of the fields at the interface between the waveguides and free 
space, and this yields an integral equation that has a well-behaved matrix 
solution. A time variation of e lujt is assumed, and all field quantities 
are in phasor form. The element locations for an x-y plane grid are 



Figc 3-2 Dielectric Covered Array of Circular Waveguides 
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given by 

-> -> 

p =mba 1 +nda 0 
mn 1 2 

as shown in Fig. 3-3 where 

b = interelement spacing along the x-direction (BX) 
d = element spacing along a^ (D) 
a = angle between a^ and a^ = grid angle (ANGLE) 
a^ja^" unit vectors 


(3-1) 


Initially an infinite array of dielectric filled circular wave- 
guide interfaced directly with free space will be analyzed. Once these 
results are obtained the extensions to allow for dielectric discs in 
the guides and dielectric cover sheets are straight forward. 


I. Waveguide fields (interior region z <_ 0) 

In this region the solutions can be expressed as the sum of a 
complete orthonormal sequence of vector modes e J } . Only the 

transverse field components are considered since when boundary conditions 
are later applied only the tangential fields are retained. Within the 
guide the transverse electric field can be expanded 

+ ^ -If. z + Y.z ^ y . z_^ 

E (int.) = Y~ (I.e J + R.e J ) %. + Y” R.e J %. (3-2) 

T J J JJ 

j=l J+l 

over the aperture 

= 0 over the rest of the cell 


where { $. } account for all TE and TM 

j mn mn 

th 

If the propagation constant of the j mode 


modes 



(vertical and horizontal) . 

a j + p ure ^y 
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imaginary (y = iB^), the mode propagates and the guide admittances are 
given by 


TE 



= Y 


o 



(3-3) 


Y 


TM 


8 . 

J 


Y 

o 



Y 

o 


(k /k ,) 
c a 


icac 

Y 


(3-4) 


where Y q = /^ = ^ and k d = - 

f d o 


The transverse H-field can be shown to be [20] 


H T (int.) 


a x y (Y I e 
z L — m m 

m 


-ik z +ik z 

- Y R e ) 0 

mm m 


(3-5) 


2. Exterior Fields (z >_ 0) 

For the region above the waveguides Floquet type modes [20,28] can 

be used to describe the fields, and for a doubly periodic planar grid 

the resultant sequence of expansion vectors are called Block functions 

± iB 2 

e 11111 } e Expanding the electric field over a grid cell in terms 

mnp 

of these vectors yields 

\ 2 

l(e xt .) -21 L L V V (3-6) 

p=l m n 

where p = 1 denotes the TE modes 
p = 2 denotes the TM modes. 
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The H-field is accordingly 

2 

H(ext.) = a x V Y Y v Y' V ( 

z l_ [_ mnp mnp mnp 

p=l m n 

For the above expansion vectors the mode admittances are given by 


Y f = Y ? 
mnl mn 


(TE) 


B BY 

mn mn o 

cay k. 


(3-8.) 


Y T = Y 
mn2 mn 


(TM) 



mn mn 


(3-9.) 


where 


Y = — and k , = ca/ye , 

o y y d 

More details on the waveguide and free space expansion functions are 
contained in subsequent sections of this chapter. 

At the aperture interface between the guide region and free space 
(z = 0) , boundary conditions require that tangential E and H be con- 
tinuous across the boundary. Hence 


E 


tangential 


= E T (int.) = E T (ext.) over the grid cell 


(3-10) 


H . - = H m (int.) = H_(ext„) over the aperture only 

tangential T T 


(3-11) 


The coefficients in the series expansion for E^ can be obtained 


using the orthogonal properties of the mode functions 
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I . 
J 


+ R 


<E(int. ) , 


^ - a h 

aperture 


d> . dA 
J 


D 

mnp 


<E (ext, ) , 


Y* > 
mnp 


« // e t 

cell 


dA 

mnp 


(3-12) 


(3-13) 


where * denotes complex conjugate. The boundary condition equation 
for tangential H is then manipulated to get an integral equation for 
tangential E field . From Eqs„ (3-5) and (3-7) 


H T (int„) 


a z x ^ (I_s - R^) Y^ <f> s 


j J j J 


(3-14) 


- H^fexto) 


= a x 
z 


III 

p m n 


D Y* * 
mnp mnp mnp 


-f 

= a x 
z 


- -tk 


y y y y ? ? n e t r da 

/ / / mnp mnp J J T mnp 

A 


p m n 


(3-15) 


From (3-15) 


V (1. - R )Y , | = V y y Y* ? // E " ? * da (3-16) 

Z_ j j j j L L L mn P mn P ' ' T mn P 


p m n 


Adding 


I 


-» -T -* 

Y 0 if E 9 <p u da 

j J J 

A 


■v 


m/ 


to both sides of (3-16) yields 


I. + R, 
j J 
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I 


I,Y.$. 
J J] 


Z. ifK ■ da 


j-i 


+ y y y y ' $ //£.$ * ^ 

i— i— L- mnp mnp J > T mn 


(3-17) 


p m n 


where Eq . (3-17) is an integral equation whose solution or unknown is the 
tangential E-field, E^* 

Substituting (3-1) for E in (3-17) gives 


i ^ y j h - i 


(I. + R.)Y. *. 
3 J J 3 


3 = 1 


+ y y y y' r // y(i. + r. )t ,'V* da 

Z_ L- l— mnp mnp i— k Tc k mnp 


p m n 


Ap . 


(3-18) 


Taking moments with respect to the set{$_,} [19] 


2 Y.I. = Y . (I . + R.) 
3 3 3 3 3 


+ 7 7 7 Y ' // ^ da // I(I+H.) T. .f * da 

i— t— mnp mnp j ^ k k k mnp 

p m n Ap. (3-19) 

The above expressed in matrix form is 


\h 





\ + r; 






I 2 + »2 

• 

Y.I. 
3 3 

= 


A 


• 

I J + R J 

0 





R J + 1 

0 





1 

1 


2 


(3-20) 


where 


A = (a„.) matrix with elements a.. 
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a.. = Y.6 . . + E E E Y' C . C * . 
ji jji pmn "“P "“PJ mnP 1 


(3-21) 


C . = // da = mode coupling coefficient 

mnpj 7 J 3 mnp 

Ap. (interior to exterior) (3-22) 


I_. = input = complex component of $ guide mode 

incident at the aperture 

->• 

= reflected component of mode 

til 

Y. = characteristic admittance of j waveguide mode 

J 

Y = characteristic admittance of radiation modes . 
mnp 


Inverting the above matrix equation gives 


I 1 + H i 


I. + R. 

3 3 


h 


= A 


-1 


2Y I 
11 


2y.I. 

J J 


(3-23) 


where 


A = inverse of A = (a 

mn 


a = element of the inverted A matrix 
mn 


The inverse can be determined on a digital computer provided the elements, 

a , of the A matrix can be found and the matrix is nonsingular. A 
mn 

discussion of convergence and truncation effects in mode coupling solutions 
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of this nature can be found in Lee et al. [24] and Mittra [25] . In this 

solution K and L which determine the number of free space modes 
max max 

were allowed to vary together and independently from 7 to 29 . No solution 

oscillation or divergence problems were observed, and since the numerical 

results also agreed with the Waveguide simulator measurements of Amitay 

and Galindo [22], accuracy was judged to be good. For many problems 

K = L =13 was used, 
max max 

3, Coupling Coefficients 

The key to practical realizations of this problem on a computer is 

that the coupling coefficients, C . , can be evaluated analytically. 
r mnpj 

As a consequence, it is unnecessary to rely on time consuming (integration 
time even on a high speed machine such as a Univac 1108 would have been 
considerable), expensive, and accuracy limiting numerical integration 
techniques. In addition, having analytic formulations of these waveguide 
to free space mode coupling coefficients, it is possible to better consider 
the number of free space modes needed to achieve reasonable convergence of 
the series for each of the A matrix elements. 

For circular waveguide apertures 


C 


rspj 


/ / 


<fr. 

3 


aperture 


a 2 tt 

Y da = / / <►. .¥ rdcj)dr 

rsp j q j rsp 


(3-24) 


where 

p = 1 denotes TE waves 
p = 2 denotes TM waves. 

The waveguide modes }are actually more complex in that (j) depends on 
the mode being TM or TE, vertical or horizontal, and marked by a two- 
dimensional wavenumber. Hence, in evaluating the coupling coefficients, 
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it is helpful to utilize a subscript and superscript notation 


c p = c p q 
rs,j rs ,mn 


Substituting for and ¥ in (3-24) and performing the weighty task 
J rsp 

of integration, it can be shown 


$ (TE) , Y (TE) 
mn rs v 


H 


2ttA i 
mn 


m -1 


k (r,s) - A 
r mn 


7 a J (X a)J (k a)/ 
2 m mn m r \ 


cos m \p 


(■ 


sin m ip 


V 


C 


* (TM), V (TM) 

mn rs 


-> 


„ . .m-l sin m ip 

2^ ^ 1 ft / 

— x 22 (-k )J (X a)J (k a) / (3 

, z , Z r m mn m r \ 
k - A \ , 

r mn cos m ip 


V V 

\ V (TE) ’ \s ( ™>> 

H N 


^ i " 1 ” 1 j (x a ) J (k a) 

k m mn m r 


( 


sin m \p 


-cos m ip 


(3 




0 (TM) , V (TE) > = 0 

mn rs' x 


H 


(3 


where 


-1 k 

V = PSI = tan 1 t*- . 

k 

x 

V denotes vertical 

H denotes horizontal. 


) 

Converting the above to the subscript and superscript notation 


(3 


c 1 2 =0 

rs ,mn 


(2 


: 1 1 
''rs ,mn 


0 . 2 . m- 1 

2tt A i 

mn a 

,2 7 N . 2 ^m v 'mn~ /w m vtv r c 

k (r,s) - A 
r mn 


J (A a) J (k a) 


r cos m \p y 
sin m ip 


V 


H 


a 


-25) 


-26) 


-27) 


-28) 


-29) 


-30) 


-31) 
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; 2 1 

^rs ,mn 


27rm i 


m-1 


J (X a)J (k a) 
m inn m r 


sin m 


-cos m ip H 


(3-32) 


0 . -2tt X' i m \ , , 

„2 2 inn r T / ^ \t/t \ 

: = o o — J ( A a ) J ( k a ) 

rs , mn .2 T 2 m mn mr 

k. — A 
r mn 


sin m \p 


cos m ip 


(3-33) 


H 


Again, the ability to evaluate these coupling coefficients without 

numerically integrating makes this an extremely effective method. In 

our computer simulation twenty-two waveguide modes were used, and the 

resultant mode index (j ) and eigenvalue for each type of mode are listed 

in Table I below. Note that the index (j) lumps information as to whether 

the mode is TE or TM and vertical or horizontal into a single index, 
mn mn 


Table I Waveguide Mode Ordering 
j $ 0 mode eigenvalues 


1 

TE n (V) 

1.87 

2 

TE U (H) 


3 

TM U (V) 

3.83 

4 

TM n (H) 


5. 

TE 21 (V) 

3.05 

6 

TE 21 (H) 


7 

tm 21 (v) 

5.13 

8 

tm 21 (h) 


9 

te 31 (v) 

4.2 

10 

TE 31 (H) 


11 

TM 31 (V) 

6.38 

12 

TM 31 (H) 


13 

te 12 (v) 

5.33 
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Table I (Continued) 

j $ 4 mode eigenvalues 


14 

TE 12 (H) 


15 

™ 12 (V) 

7.015 

16 

TM 12 (H) 


17 

TE 0 i(V) 

3 o 83 

18 

te 41 (H) 

5.31 

19 

T«01® 

2.4 

20 

™02 <H) 

5.52 

21 

te 41 (v) 

5.31 

22 

te 02 (v) 

7.015 


For most of the array configurations analyzed using this program simula- 
tion, the waveguide dimensions were chosen so that only the modes 
and $2 propagate, and hence twenty-two modes should be adequate since 
twenty of them are below cutoff. 

B. Dielectric Cover Sheets and Matching Plugs 

As has already been stated, dielectric discs within the waveguide 

can be used to generate modes or to impedance match between the guide 

fields and free space. For this purpose dielectric plugs in the ends 

of the waveguides will be considered. The program modification is 

simple in that it is only necessary to replace Y. in (14) by Y, , the 

3 3 

effective waveguide admittance at the z = 0 plane looking in the -z 
direction. Using waveguide and transmission line concepts Y. is seen 


to be 
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Y J y plug^ y 


y oj (j) + i y , (J) tan[k .. (J) * T2] 

guide plug plug 


pl„g (J) + 1 y guide (J) * tanlk plug< J) * 121 


(3-34) 


where 


v , (J) = characteristic admittance of the <f> mode in the plugs 

7 plug J 

v . , (i) = characteristic admittance of the J*-* 1 mode in the 
-’guide J 

guide (allows for dielectric fill, EPS) 

k „ (J) = phase constant of the mode in the plugs = -iy 

plug 

T2 = plug thickness , 


Two dielectric cover sheets are included in the program and should 
provide a fairly flexible basis for extending the impedance matched 
scan range of an array. The behavior of the fields in dielectric layers 
in the region z > 0 is at first not as obvious as that for the waveguide 
region- However, the fields in each of the layers and free space can be 
decomposed into TE and TM waves with respect to z [26], and then normal 
traveling wave techniques can be used to obtain the modifications of (3—19) 
needed to solve the problem. 

First, recall for systems with both electric and magnetic sources 
that the total fields are given by 

E = E A + E F = -jwA - V(V . A) - i VxF (3-35) 

k. 

H = H. + H_ = - jwF - ^ V(V • F) + - VxA (3-36) 

A F k 2 y 

Waves propagating TM and TF, to z can be determined from 

A = a (3-37) 

z 

F = a V f 
z 


(3-38) 
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where 


i p = magnetic potential 


ip = electric potential. 


For this case Eqs . (3-35)and ( 3-36 )s imp lify to [26,27] 


E = - V x(a + -r^~ V x R 

z looe 


V X (a z ^ a )j 


H = V x(a ip*) + — V x[v X (a \p f ) 1 

z loop L z J 


(3-39) 

(3-40) 


The potential functions are solutions of the homogeneous phasor wave 

equation (Helmholtz equation) subject to boundary conditions. Above each 

cell in an infinite array of the type being considered the application 

of boundary conditions gives separation constants k and k which are 

x y 

determined by the cell geometry (size and shape) and are not affected by 
the presence of the dielectric layers. Only k^ depends on the dielectric 
constant of a layer , 


\ 

/ 



Fig 0 3-4 c Array Covered with Two Dielectric Sheets. 
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The solution to the Helmholtz equation for the free space region 
(z > T3 + T4) is 


f 1 


( 

. a 


D 1 

♦ 

- z z 

pq 

,f 


D 

♦ 

p q 




z 

) 


° 


-i [k (p)x + k (p,q)y + k (p,q)z] 
x y ^ 

e 


(3-41) 


where 


k x ( p) = 


k 

o 


sin 0 cos cj) 
o o 



k (p,q) = k sin 0 sin <(> + ~ , ^P — 

y o o o d sin a b tan a 

2 2 2 2 
k Z (p,q) = k -k -k 
z v ^ o x y • 


The separation phase 


constants k and k are unchanged in a dielectric 
x y 


layer, however, 


k becomes 
z 


k z m (p,q) 



(3-42) 


e = relative dielectric constant of the m layer, 
m 

Thus it is seen that the system for z >_ 0 can be treated as a cas- 
caded waveguide- like system each section of which is characterized by a 

wave number, k . 

z 
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Fig, 3-5 . Array Covered with a Single Dielectric. Sheet. 

For a single layer of dielectric covering the array surface, the 
admittance at the z = 0 plane is 


y Q (p,q) + i y 3 (p,q) tan[k z3 (p,q) * T3] 
Y +z (p,q) “ y 3^ p,q) y 3 (p,q) + i y Q (p,q) tan[k z3 (p,q) * T3] 


(3-43) 


where 

Yo(p>q) = characteristic admittance of the ijj mode in the 
3 pq 

dielectric layer (0 z <_ T3) 
y (p,q) = free space admittance of the pqr mode (z>T3) 

k ~(p,q) = wave number of the pq^ mode in the dielectric 
z j 

T3 = dielectric thickness 

EPS3 = relative dielectric constant of the cover sheet. 

2. Two Cover Layers (EPS3, T3, EPS 4, T4) 

In this instance, the admittance looking in the +z direction at the 
level z = T3 is 


y t3 = Y T 3 (p»q> r ) = y 4 (p>q>r) 


y o (p»q»r) + i y 4 (p,q,r) tan[k z4 * T4] 
y 4 (p,q,r) + i y Q (p,q,r) tan[k z4 * T4] 


(3-44) 
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Using the above as the load for the first layer the required 
admittance at z = 0 is 


Y +Z (p> q > r ) ■ y 3 (p>q»r) 


Y t3 + i y 3 (p,q,r) tan[k z3 * T3] 
y 3 (p »q >r) + i Y t3 tan[k z3 * T3] 


(3-46) 


In the above equations 


y 0 (p>q»D 


TE 

y Q (p»q) = 


k zo (p > q) 
a) y 


k zo (p ’ q) Y o 


(3-47) 


y 0 (p»q» 2 ) 


TM coe y k 

y ° (P,q) = k zo (p ’ q) = k zo (p>q) 


k 


2 


zo 


(p»q) 




(3-48) 



[! -(T x - 




2 

2 ^q + 2 it p.. . ) i 

kd sin a kb tan or 

(3-49) 


with k being replaced by k = o)/yT = k^/y^ed to convert to the forms 
valid for the dielectrics e 


Co Waveguide Vector Modes (a = waveguide radius) 

For reference purposes this and the next section contain the 
equations defining the modes used in expanding the waveguide and free- 
space fields. In a circular waveguide of radius "a" the TE and TM 
(to z) vector modes are given by 


| (TE) = a - J (b r) 
mn r r m mn 


H 


^ . sin m<f> , 
-cos m(f) - 


+ a b J (b r) 
<p mn m mn 



(3-50) 
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V 


t mn <™> 

mn 


a (-a )J (a r) 
r mn m mn 


H 


sin mcf) 

( ) 

cos m<|> 


. ^ m T / \ 

+ a — J (a r) 
r m mn 


-cos m<|> 

( ) (3-51) 

sin m<j) 


where V denotes the vertical modes (parallel to the y-axis in Fig, 3-3) 
and H designates the horizontal modes (parallel to the x-axis). Recall 
that TE modes correspond to roots, , of 


d J (b a) 
m mn 

dr 


0 ■+ roots of J (X ) = 0 
m mn 


where X = b a. 

mn mn 

For TM modes: J (a a) = 0 -* roots X = a a 

m mn mn mn 


*Note TE (H) = 0 
on 


TH (V) = 0 
on 


D. Free Space Wave Functions (Block Functions) 

Working with the component of the wave propagating normal to the 
boundary (z-direction) 


*t y x -* 

^ (TE) k x k y 
r r 


+ik x 4- ik y 
x y 


(3-52) 


k k 

■f x -► , y -*• 

ill /m>/\ = 1 3 + 3 e 

^ (TM) k r x k r y 


ik x + ik y 
x y 


(3-53) 


2 2 2 
k -k +k 
r x y 


th 

For the mn mode and a beam pointing direction (0 ,<j> ) the 

o o 

direction cosines are 


(3-54) 
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and 


x 


T = sin 0 cos cj) 
x o T o 

(3-55) 

T = sin 0 sin $ 
y o o 

(3-56) 

\ i m 2 Tim 

= k (m) = k T — 

x ox b 

(3-57) 


k (m,n) = k T - [ — — ™ — - ■ ■ ^ | TT - — ] 

y o y d sm a b tan a 


(3-58) 


a = grid angle 

b = spacing along a^ - direction 
d = spacing along a^ direction 


grid points p = r b a- + s d a„ 


rs 


To point the beam in the (0 ,(j) ) direction, the differential phase 

o o 

shifts between elements are 


2 7T 2 7T 

6 = b T = ~ b sin 0 cos <J> 

x A x X o o 


(3-59) 


6 = d sin a T = d sin 0 sin a sin cf> (3-60) 

y X y X o o 

For z >_ 0 only the forward propagating modes are present (assumes 

system infinite in which case no reflections occur) and the Floquet 
-i 3 z 


modes are { ¥ e 
mnp 


mn 


} 


where 


mn 


= / k 2 - k 2 


o r 


o o k 2 k 2 

6 mn = k n U " ~ (~f) 1 


2 2 
= ri - f T _ 2 tutu _ /"p ^ Trn l ^ 77111 ) I 

x bx yd sin(a)k bk tan or 


(3-61) 
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Conclusions 

Grating lobe predictions are the most commonly used measure for 
estimating array scan limits, and as such the WAIM attempts in this 
study are compared with their angular predictions . As expected, it 
was observed that mutual coupling in physical arrays is such that 
reflection coefficient peaks can occur interior to the grating lobe 
angles — a condition first noted by Farrell and Kuhn [39] and a serious 
problem in any superficial design. Upon applying WAIM techniques, the 
reflection coefficient peaks have been moved outside the grating lobe 
angles thereby extending the usable scan range for a given array 
geometry. The fact that the grating lobe angle depends on the scan 
plane chosen is omitted in many discussions of grating lobes, however 
generalized grating lobe conditions will be developed here because of 
their need for some of the scan planes investigated. The grating 
lobe condition for a rectangular grid is [27] 


d = 1 

X 1 + | sin 0 m | 


(3-62) 


which can be rearranged to give 

sin 6 = -7^-1. (3-63) 

m d 

m 

For equilateral triangular grids the above equation is modified to 


d 

m 

X 


1.155 

1 +|sin 0 I* 
• m 1 


(3-64) 


Basically grating lobes appear when the grid points are phased to 
produce secondary major lobes in visible space. An alternate way of 
viewing this phenomenon is in terms of the invisible-visible mode 
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propagation approach. From this point of view, it is seen that 
grating lobes occur when a secondary free space mode makes the transi- 
tion from the invisible (nonpropagating or evanescent) to visible 
region (propagating) . The direction cosines of the pointing angle 
(6,<f>) are 

T = sin 0 cos $ (3-65) 

X 

= sin 0 sin <j> (3-66) 

For the free space region exterior to the guides 


k = k T - — ^ 
x ox b 


k = k T _ 2l ™ + 2m 

y o y d sin a b tan a 

From the separation equation, the wave numbers for the various free 
space modes are 


k (m,n) = /k 2 -k 2 (m,n) - k 2 (m,n) 
2 x y 


where 


k z (0,0) is the desired radiating mode 
2 2 2 

k > k^ + k^ -* propagating modes 

2 2 2 

k < k + k -* evanescent modes, 
x y 

th 

Since grating lobes occur when modes other than the 0 (j start to 
propagate, the critical conditions for propagation define the invis- 
ible-visible region boundaries . In direction cosine space the 
boundaries are circles such as shown in Fig. 3-5 and defined by 


k 2 = k 2 (m,n)+k 2 (m,n) = (k T - 
x y ox 


2 m* 


(k T - 
o y 


2irn 


2 m 


d sin a b tan a 


■) 


(3-67) 
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Normalized to k (3-67) becomes 
o 


€ u = (T -Si) 

r r x b 


+ (T — 

y d sin a 


m A 

b tan a 


(3-68) 


Grating lobe regions for the grid in Fig, 3-5 are the areas of mutual 
intersection which are shown shaded,, 

T 

y 


X 



Fig. 3-5 Direction Cosine Space Showing Grating Lobe Regions 

(60° grid angle) 

The following calculations show the determination of grating 
lobes angles for various scan planes for an equilateral triangular 
grid array with b = d = 1.0 and X = 1.4o The incident wave was 
polarized parallel to the y-axis in the grid plane. The spacings 
were purposefully chosen to produce a grating lobe in visible space 
(6 = 38°) so that matching effects could be evaluated, and many of 
the test programs were run using the 60° scan plane where 
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cross-polarization effects were found to be significant. The choice of 
grating lobe angle at 0 = 38° was made in an attempt to somewhat reduce 
computer time since a scan just beyond the grating lobe angle in a 
particular plane is all that is usually of interest. Thus by elimi- 
nating the necessity of scans to 90° fewer program iterations are used 
with the resultant savings of computer time. The results of all these 
simulations are contained in the Appendix. 

E-Plane Scan (9 = 90°) 


k = k , d = b = 1.0, X = 1.4, Angle = 60° 

o 


For m = 0 and n = 1, Eq. (3-68) becomes 



sin 0 


_ A l- 1 = -1 + — ■ ■ 

gl d sin a — .8666 


615 


0 


gl 


sin 1 (.615) * 38° 


H-Plane Scan (9 = 0° ) 


Here the choice n = 1 and n = 0 is made giving 


(sin - £-) + ( 


) 

b tan a 


sin 0 - lc 4 = + 



0 = 54° 


60° Scan Plane (0 = 60°) m = 1, n = 0 


(1/2 sin 0 - 1.4)2 + ( sin 0 + 
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Solving gives 0 ^ = 54° as one would expect since for the hexagonal 
grid the grid axes could be rotated 60° without changing its structure. 


30° Scan Plane (0 = 30°) 


From a symmetry consideration similar to the 0 = 60° example 


6 . = 38 c 
gl 


1. Array Performance Observations 

In the following 1^ is the applied or incident value of the 
th 

excitation of the k waveguide mode, and is the reflected component 

th th 

of the k mode. The complex reflection coefficient for the k mode 

is T k = R^/I^. Linear and circular polarizations were obtained by 
using crossed TE^ modes which correspond to modes and For 

most of these examples and are the only propagating modes in 
the waveguides, and hence are of prime interest and will be accordingly 
emphasized in plots of scan results 0 Based on many computer experi- 
mental simulations [see the Appendix for graphical results], it was 
observed that 

(1) E-plane scans were characterized by the early onset of 
grating lobe conditions (6 ^ = 38°) and probably showed the most 
intramode reflection coupling 1^ to R^. 

(2) H-plane scans yielded a grating lobe at 0 = 54° (farthest 
from broadside for this geometry) <, Match was fairly good for most 
configurations out to the grating lobe angle* Generally this was an 
uninteresting choice of scan plane. This observation should be 
carefully weighed when considering the use of waveguide simulators 
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since most simulators are compatible with H-plane scans and thus 
contribute very little scan information. 

(3) The 60° scan plane was characterized by a significant 
reflected value of the cross-polarized dominant waveguide mode 
The grating lobe was at 54° and the depolarization was found to be 
greatest near and usually just prior to the grating lobe angle,, This 
scan plane was judged to be the severest test of wide angle matching 
as a consequence of the cross-polarization effects. 

(4) The 30° scan plane had a grating lobe at 38°, but it was 
subject to less depolarization than the 60° cut. This is still a 
potentially good scan plane for initial testing. 

(5) It was difficult to get reflection peaks interior to the 
grating lobe angle with an unmodified array of circular waveguide 
elements in an equilateral triangular grid. The triangular grid and 
circular elements seem to keep mutual coupling effects to a minimum [34] 
indicating that this is a seemingly desirable type of element shape 

and array geometry 0 

(6) A high degree of coupling of power to higher order evanescent 
waveguide modes possible. In fact neglect of these modes could lead 

to some wrong conclusions. A previous study of circular aperture 
arrays [23] wrongly ignored these effects and also chose to ignore the 
problem of broadside matching which when the array is actually scanned 
lead to major design difficulties. 

(7) In choosing initial values of dielectric constants and 
thicknesses the author has proposed some approximate design methods 
which help in the optimization process [27]. 
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(8) Even though mutual coupling was low for small diameter 
elements, they were very difficult to match*, As a consequence of 
this difficulty and because of cost and production tolerance con- 
siderations at microwave frequencies, larger elements are recommended. 
The Rosenbrock optimization program was applied with waveguide 
diameter and dielectric constant as free variables, and the program 
choice for a 60° triangular grid with elements spaced d = 1.0 cm 
apart was a waveguide diameter of 0.91 cm. 

(9) Dielectric cover sheets can move reflection coefficient 
peaks interior to the grating lobe position unless properly chosen. 

2 . Computer-aided Optimization Results 

The previous observations were based on the results of many 
computer simulations selected to get physical insight into mutual 
coupling and impedance matching problems associated with an array of 
circular waveguide elements. The graphical results of these tests 
are contained Figs, A- 1 to A- 6 in the Appendix. From these runs it 
became apparent that the process of matching at broadside was diffi- 
cult enough, but to choose array parameters to increase usable scan 
range while maintaining impedance match was next to impossible by 
human trial-and-error methods. Fortunately remedies exist in the 
form of computer optimization methods, a good survey of which is 
given by Bandler [32]. The methods are generally separated into two 
groups: (1) gradient methods and (2) pattern search methods. Since 

this problem is too complex to hope gradient information could be 
determined, a pattern search technique originated by Rosenbrock[32,28] 


was chosen. 
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In brief Rosenbrock’ s method is a coordinate rotation system which 

for N variable parameters proceeds after each set of N coordinate 

searches to rotate the coordinates so that the initial increments in 

the next set of searches are along the vector direction determined by 

the previous N searches. After each successful move the step size is 

increased by a factor of three (qualitative choice of Rosenbrock) , 

and a quadratic fit is applied after the first unsuccessful move 0 The 

result is an extremely efficient pattern search in that the number of 

i 1 n 

function evaluations from an initial set of user-chosen parameter 
values is kept fairly low 0 As with all computer optimization tech- 
niques this one suffers from the fact that it cannot distinguish a 
local minimum from an absolute one. 

In this application optimizations were made first varying two 
array parameters and then varying four array parameters, A 60° scan 
plane for a hexagonal array of elements spaced l o 0 cm. apart was used 
for investigation. For the unaltered system with guides filled with 
air and having a radius of 0,44 cm, the reflection coefficient varied 
from about 0,43 at broadside to a peak of 0 o 61 at the grating lobe 
angle (0 = 54°) . The cross-polarized reflected mode showed a signifi- 
cant peak of about 0 o 3 at 0 = 48° . The main results of some of the 
optimizations are summarized as follows: 

(1) Varying the dielectric plug (T2 and EPS2) subject to a test 
2 2 

function of | R | -+- | | at 0 = 10° and 45° the Rosenbrock selections 

reduced the reflection coefficient peak (R^) to 0.37 and kept the 
reflected power in this mode under 4% over the rest of the scan range 
up to the vicinity of the 54° grating lobe where the peak was located. 


49 


(2 ) Not using matching plugs or dielectric cover sheets but 
varying guide radius (AX) - and dielectric fill (EPS) with the above 
test function, the Rosenbrock subroutine chose AX = 0.4555 and 
EPS = 1.495, With these parameters the reflection coefficient peak 
was reduced to 0.34 and reflected power was below 5.3% elsewhere. 

(3) A really significant improvement in match over the major 
portion of the scan range was obtained when four parameters were 
optimized by the Rosenbrock subroutine. In this case the dielectric 
cover sheet and the guide discs were both used with the resulting 
choices being T2 = 0.742 cm, EPS2 = 1.19, T3 = 0.127 cm, and EPS3 = 1.57. 
The test function was |r^(10°)| + |r^( 54°)|. The scan for this 
design showed a slight resonant peak of 11% reflected power at 
0 = 43°, but reflections were lowered to an almost negligible level 
of under 0.3% over the scan range away from the peak and were under 
0.3% at the grating lobe angle of 54° . 


CHAPTER IV 


EXPERIMENTAL RESULTS 


This chapter summarizes some of the experimental work done which 
has bearing on phased array matching 0 

Ac Waveguide Simulator 

There are many possible approaches to determining the performance 
of elements in a large phased array.. Since a large number of elements 
is involved, it is desirable to know the performance of various proto- 
type antennas before constructing a whole array of them. If a small 
array were being designed, measurements could be made on the array 
after construction and modifications on the individual elements might 
be achieved with little difficulty » But the number of elements in a 
large array soon makes this simple approach impractical 0 There are 
two basic techniques for modeling a large phased arrays a computer 
analysis, and simulation in a waveguide [18] > discussed here 0 

Recall from the geometry of planar phased arrays that the scan 
angle is defined as the angle of the beam maximum from broadside and 
that there are both E and H plane scans 0 The well known concept of 
component waves in a rectangular waveguide, leads to a method of simu- 
lating an infinite planar phased array o Consider the rectangular 
waveguide shown on the following page, and for simplicity, consider 
only the dominant TE^ mode whose field equations are 
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H = A cos (— x) e ^ Z 
z a 

H = i =f A sin(=) e' le * (4-1) 

X iS. ct 

c 

K 00 
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x 




Fig, 4-lo Rectangular Guide with Component Waves 0 


The expression for may be written as 


H 

z 


-2 U 


o /'^ X .oN 

-J (— + 6z) 

+ e a } 


(4-2) 


indicating the two traveling waves shown in Fig, 4-1 » Recall the 
separation equation 


k 2 =k 2 +k 2 (4-3) 

c x y 


For the case of the TE^q mode 




2 


k 


y 


= o . 


(4-4) 
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The propagation constant, 6, is given by 

2 2 2 

B = k - k Z (4-5) 

o c 

where k = cuV / y ~ 0 Thus 
o o o 


2 2 2 

k = 8 + k 

o c 


(4-6) 


Referring to Fig 0 4-1, 


From Eq e (4-7) , 


k = k sin 0 
c o 


6 = k cos 0 

o 


(4-7) 


sin 0 = 


(4-8) 


The general expression for k^ is 2n/\^ 9 where X ^ is the cutoff wave- 
length. Eq 0 (4-8) now becomes 


sin 0 = —• , (4-9) 

c 

This phenomenon in the rectangular waveguide is similar to the 
operation of a phased array elemento The angle 0 in Fig* 4-1 may be 
considered as the f, scan angle" from broadside, and various E or H 
plane scans may be simulated by exciting proper modes in specially 
built guideso For TE modes, H plane scans are simulated, while E 
plane scans are achieved by exciting TM modes <, By varying the frequency 
and guide dimensions according to Eq 0 (4-9), a number of possible scan 
conditions may be modeled „ 

The phased array simulation documented in this report was for 
a square grid of circular waveguides opening onto a conducting ground 
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plane. The geometry of the simulated array is shown in Fig. 4-2. 








Fig. 4-2 0 Geometry of Simulated Phased Array. 

Notice the waveguide cross section drawn around the central element 0 
By constructing a circular guide to be flush mounted against the 
rectangular guide and varying the frequency of the field incident on 
the circular guide, the performance of several different scan 
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conditions may be studiedo Note, however, that varying the frequency 
changes the relative dimensions of the simulated array 0 The experi- 
mental setup, illustrated in Fig c 4-3, was used to determine the 
reflection coefficient, F, for several different frequencies. 



The range of frequencies used in the experiment were chosen to 
be above the rectangular guide cutoff but below cutoff for the circular 
guide o This was done so that a matching scheme at the end of the 
circular guide was unnecessary * 

A sample calculation for a scan condition is shown below 0 


SAMPLE CALCULATION 

For a frequency of 9 GHz 

X = j = 3o 33 cm, 

For TE n mode (H-plane scan): X = 2a = 9o6 cm 

iu c 

. a X 3.33 0/c 

sin 0 = — = -T-T- = .345 

A y o O 

c 
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Therefore , 


0 = 20 o 3° (Cardinal Plane Scan) 


d 


X 


d. _ 4o 8 
X ” 3o 33 


lo 42 


£ 

X 


c 95 
3o 33 


= o 286 


Results obtained from the waveguide simulator were used to check 
the computer analysis in Chapter III, and agreement was very good. 
Table I compares the results « 


TABLE I 


Comparison 

of Computer 

and Waveguide Simulator 

Results 


Array of 

Circular Waveguides 




Reflection Coefficient 

Frequency 

Scan Angle 

Waveguide 

Computer 

(GHz) 

(H plane) 

Simulator 

Simulation 

9.0 

20.3° 

.980 

.962 

9.5 

19.2° 

.930 

.936 

10.0 

18,2° 

.890 

.900 

10.5 

17.3° 

.865 

.880 

11.0 

16.5° 

.900 

.870 

11.5 

15 o 7° 

.861 

.858 

12.0 

15.1° 

.855 

.850 


Despite the good results obtained with the simulator, it has 
several disadvantages „■ Since TE modes are most easily generated in 
rectangular guides, E plane scans are more difficult to achieve than 
H plane „ It is also difficult to get scan planes where grating lobes 


first appear. 
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Another problem, mentioned previously, is associated with the 
transition from the rectangular guide to circular guide, The diffi- 
culty arises since the cutoff frequency of the circular guide is 
usually higher than that of the rectangular guide , If the frequency 
is high enough above cutoff of the circular guide, unwanted modes 
in the rectangular guide might be excited. Operating the circular 
guide above cutoff also requires the design of a matched load 
section for terminating the circular guide. Exciting the simulator 
below circular guide cutoff presents still another problem: the 
reflection coefficients are quite large. In order to obtain meaning- 
ful results, the measurements devices must be of sufficiently wide 
dynamic range. 

Perhaps the greatest disadvantage of the waveguide simulator is 
that a particular array can be modeled at only one scan angle, A 
thorough investigation of an array would require many simulators and 
most could not be constructed with standard waveguides. The cost of 
such an analysis would be prohibitive. 

Simulation of a phased array in waveguide can be of significant 
use, however, especially the way it was employed in this report. Our 
simulator was easily fabricated, the measurements were simple, and 
the results were in excellent agreement with our computer analysis. 

In essence, our waveguide simulator was very useful, because it 
strongly supported the validity of our computer simulation, 

B, Dielectric -Cover Sheet Over a Small Array 

The impedance of a central element of a 36 element array of crossed 
slots in a ground plane was measured as a function of scan angle with 
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and without a dielectric sheet covering the array surface. The array 
operated at 1.8 GHz with crossed slots for circular polarization 
positioned as shown in Flgo 4-4, The cover sheet had a dielectric 
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Figo 4-4. Circularly Polarized Array of Crossed Slots 
(f - 1,8 GHz, d = 7 cm) 


constant of 9.0 and a thickness of 0o635 cm. The thickness was an 
arbitrarily chosen, commercially available value and was not optimally 
selected for wide angle impedance matching since such a calculation 
requires "a priori" knowledge of array impedance variations with scan 
angle. Our experimental results showed that even though introducing 
the dielectric sheet gave a slightly worse match at broadside than 
was obtained for the uncovered array, the wide angle scan matches 
were better for the dielectric covered array. 
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Co Cavity- Backed Spirals 

S-band, cavity-backed spirals were constructed and testedo Im- 
pedance versus frequency measurements showed that they had VSWR f s 
under 2 from 2„0 to 3 o 0 GHz 0 Several different baluns were built and 
tested for use with wideband spirals <, A detailed discussion of these 
baluns is found in Technical Report No 0 2 [5]« 

D. Ku-Band Spiral Array 

A four-element array of spiral radiators designed for operation 
with a center frequency of 18 o 0 GHz was photoli thcgraphically produced 
using gold on an alumina (e^ = 9) substrate c 
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Fig. A-5. E-Plane Scan for Array of Circular Waveguides with Dielectric Cover Sheet: 
a=0.44 cm, b=d=1.0 cm, X=1.4, EPS3=1.5, T3=0.1 cm, T2=0.0 
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Fig. A-6. Array of Circular Waveguides with Dielectric Cover Sheet and Dielectric Discs 
for Matching: a=-.44, b=d=1.0, X=1.4, EPS2=3.0, T2=0.2, EPS3=8.88, T3=0.02 
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Fig. A-7. Optimized Array with Cover Sheet and Dielectric Plugs for Matching 
a=0.44, b=d=l .0, EPS=1.0, EPS2=1.189, T2=0.742, EPS3=1.57, T3=0.127 
parameters chosen using Rosenbrock subroutine with TEST=R 1 (54)+R 2 (10°) 
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Fig. A-9, Array with Dielectric-Filled Waveguides: optimized (Rosenbrock) radius and 
dielectric constant, a=0.4555, EPS=1„495, b=d=l,0, x=l,4, T2=0,0, T3=0.0 
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Fig. A-10. Frequency Scanned Array of Air-Filled Guides: a = O .44 C m. , EPS = 1.0, 

b = d = 1 .0 cm. 
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Fig. A- 11. Rosenbrock Selection of Dielectric Plugs: 
X = 1.0 cm., a = 0.29 cm., b=d=0.62 cm. 
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Fig. A- 12. 60° Scan Plane: Rosenbrock Selection of Dielectric Plugs 
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